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Plasmons  describe  collective  oscillations  of  electrons.  They  have  a 
fundamental  role  in  the  dynamic  responses  of  electron  systems 
and  form  the  basis  of  research  into  optical  metamaterials1'3. 
Plasmons  of  two-dimensional  massless  electrons,  as  present  in 
graphene,  show  unusual  behaviour4'7  that  enables  new  tunable 
plasmonic  metamaterials8'10  and,  potentially,  optoelectronic 
applications  in  the  terahertz  frequency  range8'9,11,12.  Here  we 
explore  plasmon  excitations  in  engineered  graphene  micro¬ 
ribbon  arrays.  We  demonstrate  that  graphene  plasmon  resonances 
can  be  tuned  over  a  broad  terahertz  frequency  range  by  changing 
micro-ribbon  width  and  in  situ  electrostatic  doping.  The  ribbon 
width  and  carrier  doping  dependences  of  graphene  plasmon 
frequency  demonstrate  power-law  behaviour  characteristic  of 
two-dimensional  massless  Dirac  electrons4'6.  The  plasmon 
resonances  have  remarkably  large  oscillator  strengths,  resulting 
in  prominent  room-temperature  optical  absorption  peaks.  In 
comparison,  plasmon  absorption  in  a  conventional  two-dimen¬ 
sional  electron  gas  was  observed  only  at  4.2  K  (refs  13,14).  The 
results  represent  a  first  look  at  light-plasmon  coupling  in  graphene 
and  point  to  potential  graphene-based  terahertz  metamaterials. 

Graphene  is  an  attractive  two-dimensional  (2D)  carbon  material. 
Electrons  in  graphene  behave  like  massless  Dirac  fermions  and 
exhibit  many  exotic  physical  phenomena,  ranging  from  an  anomalous 
quantum  Hall  effect15,16  and  Klein  tunnelling17,18  in  electrical  transport 
to  a  universal  absorption  constant19,20  and  tunable  interband 
transitions21,22  in  optical  response.  Most  of  these  phenomena  are 
well  accounted  for  by  single-particle  excitation  of  electrons. 
Plasmons  in  graphene,  on  the  other  hand,  describe  collective 
excitations  of  2D  massless  electrons.  Previous  studies  established 
that  plasmons  have  an  important  role  in  graphene  electron  dynamics23 
and  lead  to  new  composite  excitations  in  the  form  of  plasmarons24. 
Plasmon  resonances  in  graphene  have  also  been  probed  using  inelastic 
electron  scattering  spectroscopy25,26  and  inelastic  scanning  tunnelling 
microscopy27.  However,  the  fundamental  behaviour  of  light-plasmon 
coupling  in  graphene  is  little  known,  and  the  exciting  opportunity  to 
design  terahertz  metamaterials1"3  using  patterned  graphene  structures 
has  not  been  explored. 

Here  we  report  the  first  study  of  tunable  plasmon  excitations 
and  light-plasmon  coupling  at  terahertz  frequencies  in  graphene 
micro-ribbon  arrays.  A  micro-ribbon  array  is  the  simplest  form  of 
sub-wavelength  infrared  metamaterials.  Plasmon  excitations  in  such 
arrays  correspond  to  collective  oscillations  of  electrons  across  the 
width  of  micro-ribbons.  We  show  that  such  plasmon  excitations 
can  be  controlled  through  micro-ribbon  width  (tv)  where  plasmon 
frequency  scales  as  w  1^2.  This  scaling  is  characteristic  of 


two-dimensional  electron  gases  (2DEGs)5,6,13,14  and  distinctly  different 
from  that  of  bulk  plasmons.  For  a  given  micro- ribbon  array,  plasmon 
resonance  can  be  further  tuned  in  situ  using  electrostatic  gating. 
Plasmon  frequency  varies  with  carrier  concentration  (w)  as  n1^4,  a 
signature  of  massless  Dirac  electrons4"6.  Comparisons  between  exper¬ 
imental  plasmon  frequencies  and  simulation  results  yield  an  effective 
dielectric  screening  k  of  5  in  our  devices.  The  observed  light- 
plasmon  coupling  from  massless  electrons  in  graphene  is  remarkably 
strong:  the  integrated  oscillator  strength  of  graphene  plasmons  is  an 
order  of  magnitude  larger  than  that  achieved  in  conventional  2D 
electron  systems13,14.  Hence,  graphene  plasmon  excitations  lead  to 
prominent  terahertz  absorption  peaks  at  room  temperature,  but  low 
temperatures  (4.2  K)  were  required  to  measure  plasmon  absorption 
in  conventional  2DEGs13,14.  The  strong  and  tunable  plasmon-light 
interaction,  together  with  the  excellent  electrical  properties  of 
graphene  and  its  compatibility  with  micro/nanofabrication,  holds 
great  promise  for  graphene-based  terahertz  metamaterials  and 
optoelectronic  devices. 

Electromagnetic  radiation  cannot  couple  to  plasmon  excitations 
in  homogeneous  graphene  sheets,  but  can  readily  excite  plasmon  res¬ 
onances  in  engineered  graphene  structures  with  dimensions  much 
smaller  than  light  wavelengths.  In  our  study  we  used  a  simple  form 
of  metamaterial  made  up  of  periodic  graphene  micro-ribbon  arrays 
illustrated  in  Fig.  la  (not  to  scale)  for  terahertz  plasmon  excitations. 
Ribbon  width  varied  from  1  to  4  |xm  and  gap  and  ribbon  width  ratio 
was  maintained  at  1:1.  We  controlled  the  carrier  concentration  in 
graphene  ribbon  arrays  using  an  ion-gel  top  gate28,  which  allows  a 
large  doping  range  through  electrostatic  gating.  Figure  lb  shows 
the  schematic  side  view  of  a  typical  gated  device.  The  atomic  force 
microscopy  (AFM)  image  of  a  micro-ribbon  array  sample  with  a 
4  pm  ribbon  width  and  its  gate-dependent  electrical  resistance 
trace  are  shown  in  Fig.  lc,  d,  respectively.  The  device  has  a  charge 
neutral  point  (CNP)  at  0.2  V  top  gate  voltage  (Vg),  and  our 
samples  have  mobilities  around  1,000  cm2  V  1  s  \ 

We  studied  plasmon  resonances  and  light-plasmon  coupling  in 
graphene  micro-ribbon  arrays  using  polarized  Fourier  transform 
infrared  spectroscopy.  Infrared  absorption  from  doped  charge  carriers 
is  directly  related  to  —  AT/Tcnp,  the  gate-induced  decrease  of 
transmission  through  graphene  micro-ribbon  arrays.  Here  TCNP  is 
the  transmission  coefficient  at  CNP  and  AT  =  T  —  TCNP.  The 
— AT/Tcnp  spectra  of  a  gated  4  |xm  sample  (Tg=  -2.0  V)  for 
terahertz  radiation  polarized  parallel  and  perpendicular  to  the 
ribbons  are  shown  in  Fig.  le  and  f,  respectively. 

For  incident  light  polarized  along  the  ribbons,  the  electrical 
response  of  charge  carriers  is  similar  to  that  of  free  electrons  in  a 
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Figure  1  |  Plasmon  resonance  in  gated  graphene  micro-ribbon  arrays,  a,  Top-view  illustration  of  a  typical  graphene  micro-ribbon  array.  The  array  was 
fabricated  on  transferred  large-area  CVD  graphene  using  optical  lithography  and  plasma  etching,  b,  Side  view  of  a  typical  device  incorporating  the  graphene 
micro-ribbon  array  on  a  Si/Si02  substrate.  The  carrier  concentration  in  graphene  is  controlled  using  the  ion-gel  top  gate,  c,  AFM  image  of  a  graphene 
micro-ribbon  array  sample  with  a  ribbon  width  of  4  pm  and  a  ribbon  and  gap  width  ratio  of  1:1.  d,  Gate-dependent  electrical  resistance  of  this  graphene 
micro-ribbon  array.  The  resistance  has  a  maximum  at  charge  neutral  point  \/CNP  =  0.2  V.  e,f,  Gate-induced  change  of  transmission  spectra,  —  (T  -  TCNP)/TCNP 
(red  solid  line),  with  incident  light  polarized  parallel  (e)  and  perpendicular  (f)  to  the  ribbon  length,  respectively.  The  gate  voltage  was  set  at  Vg=  —2  V.  For 
parallel  polarization  in  e,  the  response  originates  from  free  carrier  oscillation  and  can  be  well  reproduced  by  a  Drude  fit  (black  dashed  line).  For  perpendicular 
polarization  in  f,  the  spectrum  shows  a  prominent  absorption  peak  at  3  THz  (1  THz  =  33.3  cm-1)  because  of  plasmon  excitation.  Plasmon  resonance  is 
characterized  by  a  Lorentzian  lineshape  (blue  dashed  line).  A  small  free  carrier  contribution  described  by  Drude  absorption  (magenta  dashed  line)  is  also 
present  as  a  result  of  graphene  absorption  outside  the  fabricated  micro-ribbon  array  area.  The  plasmon  absorption  of  over  13%  is  remarkably  strong,  and  its 
integrated  oscillator  strength  is  more  than  an  order  of  magnitude  larger  than  that  achieved  in  2DEGs  in  conventional  semiconductors. 


homogeneous  graphene  sheet,  and  the  spectral  shape  is  well  described 
by  the  Drude  model  as  Im(—  l/(<u+  rTD))  (black  dashed  line  in 
Fig.  le)  with  a  scattering  rate  rD  of  4  THz.  This  Drude  absorption 
decreases  monotonically  with  frequency  co.  For  incident  light 
polarized  perpendicular  to  the  ribbon,  the  resulting  spectrum  is 
completely  different:  an  absorption  peak  originating  from  plasmon 
oscillation  dominates  the  optical  response  (Fig.  If).  This  plasmon 
absorption  lineshape  can  be  described  by  a  damped  oscillator  as 
Im(  —  <u/(w2  —  (t>2  T  iruTp))  (blue  dashed  line)  with  a  plasmon  reson¬ 
ance  (Up  of  3  THz  and  a  spectral  width  Fp  of  4  THz.  There  is  also  a  small 
background  contribution  from  the  free  carrier  absorption  (magenta 
dashed  line)  because  part  of  the  infrared  beam  transmits  through 
unpatterned  graphene  around  the  micro-ribbon  array.  The  light— 
plasmon  coupling  in  graphene  is  remarkably  strong,  with  over  13% 
absorption  at  the  plasmon  resonance  for  one  doped  carbon  monolayer. 
The  observed  plasmon  resonance  linewidth  is  similar  to  the  Drude 
scattering  rate  for  parallel  light  polarization,  suggesting  that  plasmon 
broadening  is  largely  limited  by  the  same  scattering  processes  affecting 
electrical  transport. 

The  light-plasmon  coupling  in  graphene  is  remarkably  strong 
compared  with  that  in  conventional  2DEGs  in  semiconductors,  and 
arises  from  both  the  small  effective  electron  mass  and  the  efficient 
electrical  gating  in  graphene.  In  conventional  2DEGs,  the  field- 
induced  carrier  concentration  is  limited  to  ~1  x  1012  cm  2  to  avoid 
semiconductor  dielectric  breakdown.  At  this  carrier  density,  the 
effective  mass  of  an  electron  at  Fermi  energy  is  0.02  m0  in  graphene. 
This  electron  effective  mass  is  3.5  (13)  times  smaller  than  that  in 
GaAs  (Si),  and  the  integrated  plasmon  oscillator  strength  is  corre¬ 
spondingly  larger.  In  addition,  graphene  can  be  gated  much  more 
efficiently  to  over  1  x  1013  cm  2  because  of  its  excellent  chemical 
stability  and  compatibility  with  different  gating  configurations. 


The  resulting  light-plasmon  coupling  in  graphene  can  be  an  order 
of  magnitude  stronger  than  that  achieved  in  conventional  2DEGs, 
and  enabled  us  to  observe  terahertz  plasmon  resonance  for  the  first 
time  at  room  temperature. 

The  plasmon  resonance  in  graphene  micro-ribbon  arrays  can 
be  tuned  in  situ  by  electrostatic  doping.  For  simplicity,  we  will 
focus  on  the  hole-doping  regime.  We  determine  the  gate-induced 
Fermi  energy  shift  and  carrier  concentration  in  graphene  through 
interband  optical  transitions21,22.  This  is  based  on  the  effect  that 
interband  transitions  can  be  blocked  up  to  2|EF|  energy  in  hole- 
doped  graphene  as  a  result  of  empty  initial  states.  Figure  2a  shows 
gate-induced  changes  of  interband  transmission  spectra  in  the 
3,000-9,000  cm  1  spectral  range  at  different  gate  biases.  Increased 
infrared  transmission  is  observed  up  to  a  threshold  energy  that 
shifts  to  higher  values  with  increased  doping.  This  threshold 
energy  is  at  2|EF|  (refs  21,22),  from  which  carrier  concentration 
can  be  deduced  by  n  =  (|£f|/^vf)2/tt  (refs  15,16).  Plasmon  exci¬ 
tations  in  the  graphene  micro-ribbon  array  can  be  varied  by 
electrical  gating,  as  shown  in  Fig.  2b  for  absorption  of  perpendicu¬ 
larly  polarized  light.  Prominent  plasmon  absorption  peaks  are 
observed,  which  shift  to  higher  energies  and  gain  oscillator  strength 
with  increased  carrier  concentration.  For  comparison,  we  show 
(in  the  inset  of  Fig.  2b)  free  carrier  absorption  spectra  at  different 
gate  voltages  (probed  with  light  polarized  parallel  to  micro¬ 
ribbons).  The  spectral  shape  remains  the  same  in  all  spectra, 
although  the  oscillator  strength  increases  with  carrier  concentration. 

Plasmon  excitations  in  graphene  micro-ribbon  arrays  can  also  be 
controlled  by  engineering  the  ribbon  width.  Figure  2c  shows  AFM 
images  of  three  micro-ribbon  arrays  with  1,  2  and  4  |xm  ribbon 
widths,  respectively,  and  Fig.  2d  shows  —  AT/TCNP  spectra  for  the 
three  structures.  All  spectra  were  taken  at  the  same  hole-doping 
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Figure  2  |  Control  of  plasmon  resonance  through  electrical  gating  and  micro-ribbon  width,  a,  Mid-infrared  transmission  spectra,  T/TCNP,  of  the  graphene 
ribbon  array  in  Fig.  1  as  gate  voltage  Vg  —  VCNP  varies  from  —0.3  to  —2.2  V.  The  voltages  corresponding  to  the  unlabelled  lines,  starting  with  the  red  line  and 
alternating  downwards,  are:  —2.0  V,  —1.6  V,  —1.2  V,  —0.7  V  and  -0.3  V.  On  electrical  gating,  optical  transmission  is  increased  up  to  a  threshold  energy 
of  2|Ef|  as  a  result  of  blocked  interband  optical  transitions.  This  threshold  energy  provides  direct  determination  of  Fermi  energy  EFand  carrier  concentration 
n  =  (\EF\/tlvF)2/'n  in  gated  graphene,  b,  Control  of  terahertz  resonance  of  plasmon  excitations  through  electrical  gating.  Terahertz  radiation  was  polarized 
perpendicular  to  the  graphene  ribbons.  The  plasmon  resonance  shifts  to  higher  energy  and  gains  oscillator  strength  with  increased  carrier  concentration.  For 
comparison,  the  inset  shows  corresponding  spectra  due  to  free  carrier  absorption  for  terahertz  radiation  polarized  parallel  to  the  ribbons.  For  this  polarization, 
absorption  strength  increases  with  carrier  concentration,  but  spectral  shape  remains  the  same,  c,  AFM  images  of  samples  with  micro-ribbon  widths  (w)  of  1, 

2  and  4  pm.  d,  Change  of  transmission  spectra  with  different  graphene  micro-ribbon  widths  for  the  same  doping  concentration  of  1.5  x  1013  crrT2  The 
Drude  background  contributed  by  unpatterned  graphene  around  the  arrays  (as  in  Fig.  If)  was  subtracted,  and  all  spectra  were  normalized  by  their  respective 
peak  values  for  convenience  of  comparison.  Plasmon  resonance  frequency  cop  shifts  from  3  to  6  THz  when  micro-ribbon  width  decreases  from  4  to  1  pm. 

concentration  of  1.5  x  1013  cm  2  and  were  normalized  by  the 
respective  peak  values  for  convenience  of  comparison.  The  plasmon 
resonance  clearly  shifts  to  higher  energy  with  decreasing  ribbon 
widths:  resonance  frequencies  are  at  3,  4.1  and  6  THz  for  ribbon 
arrays  of  4,  2  and  1  pm  ribbon  widths,  respectively.  Using  nano- 
and  micro-ribbon  arrays  with  a  larger  range  of  ribbon  widths,  tera¬ 
hertz  resonance  from  1-10  THz  can  be  straightforwardly  engineered. 

Next  we  examine  the  quantitative  scaling  behaviour  of  plasmon 
resonance  frequency  with  carrier  concentration  and  ribbon  width. 

Figure  3a  plots  plasmon  resonance  frequencies  as  a  function  of  |£F| 

(or  equivalently  |n|l,/2  in  the  top  label)  for  different  micro-ribbon 
widths.  For  a  given  carrier  concentration,  the  plasmon  resonance 
frequencies  in  different  graphene  micro-ribbon  arrays  are  described 
by  a  scaling  behaviour  of  cop  oc  w~1/2.  This  is  best  shown  in  Fig.  3b, 
where  all  data  in  Fig.  3a  fit  into  a  universal  curve  once  the  plasmon 
resonance  frequencies  are  rescaled  by  w~1/2.  This  scaling  of  cop  with 
the  inverse  square  root  of  the  relevant  spatial  dimension  is  character¬ 
istic  of  2DEGs5,6,  in  contrast  to  a  constant  plasmon  resonance 
frequency  in  bulk  plasmons.  The  universal  relation  between  cop  and 
|£p|  in  Fig.  3a  can  also  be  described  by  a  power-law  scaling  with  top  oc 
Ify1/2  oc  rc1/4  (solid  line).  This  power-law  dependence  of  n1^4  is  a  sig¬ 
nature  of  massless  Dirac  electrons4-6,  which  is  in  contrast  to  « 1/2  scaling 
(dashed  line)  in  conventional  semiconductors. 

A  theoretical  understanding  of  the  scaling  behaviour  of  plasmon 
resonance  in  micro-ribbon  arrays  can  be  achieved  using  the  quasi¬ 
static  description,  which  is  a  good  approximation  for  ribbon  widths 
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(~1  pm)  much  smaller  than  photon  wavelengths  (~100  pm).  The 
Maxwell  equations  satisfying  the  2D  ribbon  array  boundary  conditions 
can  be  simplified  to 


9 'y  ku>  dx  \  dx 


(1) 


where  k  is  the  effective  dielectric  screening  constant,  cr2D  is  the  2D 
high-frequency  conductivity  at  carrier  concentration  n,  and  F(x)  is 
a  square  function  with  F(x)  =  1  at  ribbon  positions  and  F(x)  =  0 
in  the  gap  (see  Supplementary  Information).  By  introducing 
dimensionless  parameters  X=x/w  and  Y  =  y/w  and  using  the 
expression  of  cr2D  for  massless  Dirac  electrons,  the  above  equation 
can  also  be  expressed  as 


9d>  9  /  9<F> 

dY  =  r}dx{F(x)ax 


(2) 


with  the  dimensionless  parameter 


i?  = 


T  1/2  2  1/2 

2^ t  '  vFe  n 
(o(co  +  iT)  h  kw 


The  electrodynamic  responses  of  a  micro-ribbon  array  are  uniquely 
determined  by  the  dimensionless  parameter  17.  Maximum  absorption 
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Figure  3  |  Scaling  laws  of  graphene  plasmon  resonance  frequency,  a,  Plasmon  resonance  frequency  wp  as  a  function  of  |EF|  (or  equivalently  |n|1/2  in  the  top 
label)  for  micro-ribbon  arrays  of  different  widths,  b,  Plasmon  excitation  <up  was  normalized  by  1  /Jw  for  micro-ribbon  arrays  of  different  widths,  which  fits  all 
data  points  (symbols)  into  a  universal  curve  (solid  line).  This  w-1/2  scaling  of  wp  is  characteristic  of  2D  electron  systems.  The  universal  doping  dependence 
of  plasmon  resonances  is  described  by  a  scaling  law  of  wp  oc  |EF|^2  oc  n1/4  This  n1^4  scaling  law  is  a  signature  of  massless  Dirac  fermions.  In  comparison,  wp 
scales  as  n1/2  (dashed  line)  in  conventional  semiconductors. 


a 


X  (pm) 


X  (pm) 


X  (pm) 


Figure  4  |  Simulation  of  plasmon  excitations,  a,  Transmission  change  spectrum  —AT/TCNP  simulated  by  finite  element  analysis  (dashed  line)  for  the  sample 
in  Fig.  1  at  carrier  concentration  of  1.5  x  1013  citT2  It  reproduces  well  the  experimentally  observed  spectrum  (solid  line)  when  the  effective  environment 
dielectric  constant  k  was  set  as  5,  and  corresponds  to  an  electron-electron  interaction  strength  of  e2//c/ivF«;  0.4.  b-d,  Simulation  results  for  current  density 
amplitude  (upper  panel)  and  phase  (lower  panel)  of  the  device  at  frequencies  below  resonance  (1  THz,  b),  at  resonance  (3  THz,  c)  and  above  resonance 
(5  THz,  d).  The  charge  carriers  oscillate  perpendicular  to  the  graphene  ribbons  on  terahertz  irradiation.  The  oscillating  current  is  highest  at  the  plasmon 
resonance  frequency.  The  relative  phase  of  the  oscillating  current  with  reference  to  the  incident  electrical  field  also  varies  quickly  and  changes  sign  at  the 
resonance  frequency.  Both  are  characteristics  of  a  resonant  excitation. 


of  incident  light  takes  place  at  a  specific  17  value,  which  determines 
the  plasmon  resonance  frequency  <wp  for  given  values  of  w  and  n. 
When  either  ribbon  width  or  carrier  concentration  changes, 
plasmon  frequency  will  change  accordingly  following  the  scaling 
behaviour  of  u>v<^nl^-w  *'2  (for  weak  damping  rate  T),  as 
observed  experimentally. 

Although  the  differential  equation  readily  explains  the  scaling 
behaviour  of  plasmon  resonance  frequency,  it  cannot  be  solved 
analytically  for  a  quantitative  value  of  a>p.  Instead,  we  use  numerical 
calculations  based  on  finite  element  methods  to  simulate  this 
system,  which  has  only  one  unknown  parameter,  namely  the  dielec¬ 
tric  constant  k.  Using  a  k  value  of  5,  the  simulation  (dashed  line  in 
Fig.  4a)  reproduces  well  the  observed  —  AT/TCNP  spectrum  for  a 
ribbon  width  of  4  p,m  and  a  doping  concentration  of  1.5  x 
1013  cm  2  (solid  line  in  Fig.  4a).  The  k  value  of  5  is  reasonable  in 
our  devices,  given  that  the  terahertz  dielectric  constants  of  ion  gel, 


Si02  and  silicon  are  3  (ref.  29),  4  and  11  (ref.  30),  respectively. 
Such  a  dielectric  screening  constant  will  lead  to  an  electron-electron 
interaction  strength  of  aG  =  e2/ Khvp  «  0.4. 

Figure  4b-d  plots  simulation  results  for  current  density  amplitudes 
(top  panels)  and  phases  (bottom  panels)  in  the  same  micro-ribbon 
array  for  frequencies  below  resonance  (1  THz),  at  resonance  (3  THz) 
and  above  resonance  (5  THz),  respectively.  The  current  distribution 
describes  electrons  oscillating  back  and  forth  across  the  width  of 
graphene  micro-ribbons  on  terahertz  radiation.  The  oscillating 
current  amplitude  is  largest  at  the  plasmon  resonance  frequency. 
Simultaneously,  the  phase  of  the  current  oscillation  relative  to  that  of 
the  incident  electric  field  varies  quickly  and  changes  sign  across  the 
resonance  frequency.  Both  behaviours  are  typical  characteristics  of 
resonance  excitation. 

Our  studies  show  that  plasmon  resonances  in  graphene  can  be 
controlled  in  structure-engineered  graphene  micro-ribbon  arrays. 
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which  represent  the  simplest  form  of  graphene  terahertz  metamater¬ 
ials.  Plasmons  couple  to  terahertz  radiation  strongly,  absorbing  over 
13%  at  the  plasmon  resonance  for  half  a  monolayer  material  at  room 
temperature.  Peak  absorption  at  resonance  can  become  even  stron¬ 
ger  with  reduced  plasmon  linewidth,  which  is  currently  limited  by 
the  scattering  rate  as  in  Drude  conductivity.  We  expect  that  in 
higher-quality  graphene  samples,  the  plasmon  linewidth  will  be 
reduced  in  proportion  to  the  increase  in  mobility.  The  strong  and 
broadly  tunable  plasmon  resonances  observed  here  in  graphene 
micro-ribbon  arrays  can  be  generalized  for  designing  more 
sophisticated  metamaterials  based  on  graphene,  and  enable  new 
ways  of  controlling  terahertz  and  far-infrared  radiation. 

Methods 

Large-area  graphene  was  grown  using  chemical  vapour  deposition  on  copper  and 
transferred  the  graphene  sheet  onto  a  Si/Si02  wafer  following  ref.  31.  Au/Cr 
electrodes  of  50  nm  thickness  were  deposited  onto  the  graphene  sheet  in  vacuum 
using  stencil  masks.  Micro-ribbon  arrays  were  fabricated  in  a  2.5  x  2.5  mm2  area  at 
the  centre  of  the  graphene  sheet  using  standard  optical  lithography  followed  by 
oxygen  plasma  etching. 

Infrared  transmission  measurements  were  performed  with  a  Fourier 
transform  infrared  spectrometer  at  the  Advanced  Light  Source  in  Lawrence  Berkeley 
National  Laboratory.  All  measurements  were  performed  in  0.1  torr  vacuum  at 
room  temperature. 
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